Spin superconductor in ferromagnetic graphene 
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We show a spin superconductor (SSC) in ferromagnetic graphene as the counterpart to the charge 
superconductor, in which a spin-polarized electron- hole pair plays the role of the spin 2{h/2) 'Cooper 
pair' with a neutral charge. We present a BCS-type theory for the SSC. With the 'London- type 
equations' of the super-spin-current density, we show the existence of an electric 'Meissner effect' 
against a spatial varying electric field. We further study a SSC/normal conductor/SSC junction 
and predict a spin-current Josephson effect. 
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The superconductivity was discovered about a century 
ago.[l| Sirice then, it has been one of the centrai sub- 
jects in physics. [1] Many fascinating properties of su- 
perconductors, such as zero resistance,[H the Meissner 
effect, [2] and the Josephson effect, [3] have many appli- 
cations nowadays. On the other hand, the potential ap- 
plication of the spin degrees of freedom of an electron, 
the fìeld of spintronics, is stili rapidly developing and is 
emerging as a major field in condensed matter physics. [4] 

The key physics of superconductivity was well under- 
stood in the BCS theory: [5] Electrons in a solid state 
system may have a net weak attraction so that they 
form Cooper pairs which can then condense into the BCS 
ground state. The simplest 5-wave Cooper pairs are of 
electric charge 2e and spin singlet. A dual of super- 
conductor is the so-called exciton condensate in which 
a Cooper pair-like object is a particle-hole pair which 
is charge-neutral while its spin may either be singlet or 
triplet. We name a spin-triplet exciton condensate as 
the spin superconductor (SSC). The exciton condensates 
can exist in many physical systems.[6|,0] However, a gen- 
eral drawback of the exciton condensate is its instability 
because of electron-hole (e-h) recombination that lowers 
the total energy of the system (see Fig.l(a)). The typ- 
ical lifetime of an exciton is restricted from pico-second 
to nano-second, and at most limited to the micro-second 
range.jsl-llQl It is too short for many meaningful applica- 
tions. Thus, finding a long-lived exciton gas becomes an 
important task. 

Recently, graphene, a single-layer hexagonal lattice of 
carbon atoms, has been successfully fabricated. |1 ll-ll4 1 
The unique structure of graphene leads to many pecu- 
liar properties. e. g., the relativistic-like quasi-particles 
spectrum.0, ^] For graphene, the charge carriers are 
usually spin-unpolarized. However, if graphene is grow- 
ing on a ferromagnetic (FM) material[15-17] or is under 
an external magnetic fìeld[18], a spin split M can be in- 
duced. Then the carriers are spin-polarized and the Dirac 



points with different spins are split. When the Fermi 
level lies in between the spin-resolved Dirac points (see 
Fig. HJb)), the spin-up carriers are electron-like while the 
spin-down ones are hole-like. These positive and nega- 
tive carriers attract and form e-h pairs that are stable 
against the e-h recombination due to the Coulomb inter- 
action. This is because the filled electron-like states are 
now below the hole-like states, as shown in Fig. HJ^b), 
unlike in conventional exciton systems in semiconductors 
(Fig. m^a)) where the electron states are above the hole 
states. If a carrier jumps from the electron-like state to 
the hole-like one, the total energy of the system rises. 
This prevents the e-h recombination and means the e-h 
pairs in FM graphene is stable and can exist indefìnitely 
in princeple. Therefore, this e-h pair gas can condense. 
In this work, we show that this condensate is a SSC and 
the spin current is dissipationless. We derive the London- 
type equations of the super-spin-current and fìnd an elec- 
tric 'Meissner effect' against the spatial variation of an 
electric field. 

We consider an interacting electron system in graphene 
with the Hamiltonian H = Hq -\-Uc where Hq is the free 
Dirac fermion Hamiltonian and Uc the electron-electron 
(e-e) Coulomb interaction: 
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where ^kcr = (<^kcr7 ^kcr)"^? S\^c7 = cl-,^) are the Fourier 
components of the electron annihilation operators Si^y at 
sites i for the sublattices 5, and = s\^Sia- are the 
locai electron number operators. k = (kx^ ky) is the mo- 
mentum, a = (t,i) represents the spin, M is the FM 
exchange split energy, Uf^ is the e-e Coulomb potential, 
and vf = '^ta^/^ with the nearest hopping energy t and 
the carbon-carbon distance a^. Here we have ignored the 
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valley degree of freedom, because the two valleys are de- 
generate and the inter-vahey couphng is normany very 
weak due to the two valleys being well separated in k- 
space. Hereafter we also set the Fermi energy Ep at zero. 
By taking a unitary transformation: aka = ^^TC^ai^ra 
and òkcr = cakrcr with the pseudo spin index r = ±, 
c = e*^/^/>/2 and = tan~-^(/c^//ca^), the free Hamilto- 

nian Hq can be diagonalized Hq = Z]k,r,cr ^^^^kraCi^kra' 
where e^-cr = —crM -\- rvpk are four energy bands (see 
the blue curves in Fig.ld) because the spin-degeneracy is 
Hfted now. While Ef = 0, e-^ and are high-energy 
bands. In the following, we focus on the low energy part 
and only two bands and e_| are involved. Here the 
band is electron-like, while the band is hole-like 
and the annihilation operator o^k-j, also means to create 
a spin-up hole. Thus, we defìne operators aket — <^k+t 
and a^.^^ = ak-|. The Hamiltonian Hq can then be 
written as: 

For the e-e interaction /7c we also focus on the 
two low-energy bands which are given by the terms 

^k-q,et^keto^k'+q,/it^k'^f Furthermore, we keep only 
the terms whose mo menta satisfy k — q = k', giving rise 
to the zero momentum e-h pair that is energetically fa- 
vor able. Under these approximations, the interaction Uc 
reduces to the attraction between electrons and holes 

^^C = - ^ ^kk'Q^k/et^k'/it^^^t^^e^' 
k,k' 

where /7kk' = {U^l^e'^^'-^^ + U^^e'^^'^'^ + + 
^kkO/4 with /7£^, = /7aj>e-^(^-^')-(^^+^) and U^^ = 
/7o|e~*^^~^ "^'^^ for the coordinate vj of the site j and 
the lattice spacing vector S. /7kk' is a large positive vaine 
at k = k' and it gradually and oscillatorily decays to 
zero with increase of |k — k'|. As discussed before, this 
attractive interaction does not induce the e-h recombi- 
nation while it binds the electrons and holes into pairs. 
The mean field approximation of Eq. (|3]) reads 

k k 

with the e-h pair condensation order parameter 

Ak =- X]k' ^k'k(<^k'/it<^k'et)- 

Comparing with the spin singlet Cooper pair with 
charge 2e, this e-h pair is of spin h and charge neutral. 
The total mean field Hamiltonian then is given by 




The energy spectrum for the mean field Hamiltonian 
Hmf is shown in Fig. Hl^d). An energy gap with the 



magnitude of |Ak| is opened. When an electron and a 
hole combine into an e-h pair, the energy of the system 
is reduced by 2|Ak|. This means the condensed state of 
the e-h pairs is more stable than the unpaired one. Thus, 
the ground state of FM graphene is a neutral superfiuid 
with spin h per pair, namely, a SSC state. The spin cur- 
rent can dissipationlessly flow in the SSC and its spin 
resistance is zero. 

The energy gap A can be estimated as follows. By us- 
ing the definition Ak = — Xl^' ^k'k(c»^k'^t^k'et) the 
Hamiltonian Hmf^ one has the self-consistent equation: 

Ak = ^(f/kk'AkV2^){/(-A) - /(A)}, (5) 

k' 

w here f{A) l/[exp(A/A:BT) + 1], A = 

Y^(M - k'Y + Ak' and T is the temperature. At 
zero temperature and assuming t/kk' = UO{kD — |k — k'|) 
with the cut-off momentum ko-, the self-consistent 
equation ([5j) reduces to 1 = {U/2)^^6{kD — 
k)/J(M-ky^A^ = ^ C deu , '"^^ 

where = vpkjj. Numerically, we solve the self- 
consistent equation by using the e-e interaction Uc 
with the nearest neighbor cut-off. The gaps vary as the 
cut-off Cd for a fixed M or as the FM split energy M for 
a fixed Cd are shown in Fig. Hl^e) and (f), respectively. 
We see that the gap A grows faster than an exponential 
function with increase of cd- When M = bmeV and 
eo O.18t,0,|Ì^ one gets A ^ 3meV. This yields the 
criticai temperature Te of the transition from the normal 
state to the SSC at about 30i^. Similar to the case in a 
super conductor, in the presence of weak impurities. Te 
is slightly reduced but the SSC phase can stili exists if 
T < Tc^ except in the case when the impurity strength 
is larger than A and its density is higher than 1 /^^ with 
^ being the coherence length ^ = Hvf/ A- 

Meissner effect is the criterion that a superconductor 
differs from a perfect metal: The magnetic field can not 
enter the bulk of a superconductor. [2] This phenomenon 
can be described by the London equations.[19] Is there 
a Meissner- like effect for the SSC? Consider a SSC with 
the superfiuid carrier density in an electric field E 
and a magnetic field B. A magnetic force F = (m • V)B 
acts on these spin carriers. Here m = {AiigiiB /h)s is the 
magnetic moment of a carrier; /i^ is the Bohr magneton 
and g is the Lande factor. This force accelerates the 
carrier by Newton's second law F = m*dv /dt for a carrier 
with the velocity v and the effective mass m*. The spin 
current density JJg = n^vs is thus a tensor. The time 
derivative of this super-spin-current density JJg is then 
given by 

d^s/dt = a{s ■ V)Bs, (6) 

with the Constant a = ATrg/UBTis/hm^ . Comparing 
with London's first equation for the super-charge-current 
density: [19] dJ/dt oc E, the spatial variation of B along 




FIG. 1: (color online). The schematic diagrams: (a) the band 
structure for the conventional exciton system; (b) for FM 
graphene; and (c) for the proposed four-terminal device to 
measure the SSC state, (d) The schematic energy bands of 
FM graphene for free electrons (blue curves) and for the e-h 
pair condensate (black curves). (e) The gap A vs. cd at the 
FM magnetic moment M — òmeV and (f) The gap vs. M at 
CD 0.18t. 

with the magnetic moment plays the role of an external 
field accelerating the spin carriers. 

When an electric fìeld E applies, by acting (s-V) on two 
sides of the Maxwell equationsV x B = ji^e^d^/dt and 
using Eq.®, we obtain ^[V x JJ^] = ^[ji^e^ais • V)Es]. 
Integrating over the time one has the equation for JJ^ 

V X JJ5 = /ioeo«(s • V)Es, (7) 

where the integrai Constant is taken to be zero because of 
the requirement of thermodynamic equilibrium. Instead 
of the magnetic field in London's second equation for the 
superconductor,[19] the 'external field' here is the spatial 
variation of the electric field E along with the magnetic 
moment. 

Eqs. (|6|) and ([7]) for JJ5 play roles similar to the Lon- 
don equations in super conductor.[l9j For example, if the 
system is in the steady state, dìg/dt = a(s • V)Bs = 
implies that the variation (s • V)B of the magnetic field 
along the direction of the FM magnetic moment m must 
be zero because of the zero spin resistance. On the other 
hand, Eq.([7j) means the variation of the electric field, 
(s • V)E, is zero in bulk of the SSC. This is an electric 
'Meissner effect' in the SSC against a spatial variation of 
an electric field. 

We now give an example of this electric Meissner effect. 
Consider a positive charge Q at the origin and an infinite 
FM graphene in the x-y piane at z = Z as shown in Fig. 
[2fa). The charge Q generates an electric field E in FM 
graphene piane. This electric field will induce a super- 
spin-current in graphene against the spatial variation of 
E. Assuming that the magnetic moment m (i.e., s) is in 




FIG. 2: (color online), (a) The schematic diagram for the 
device consisting of a positive charge Q and FM graphene. (b) 
The variation dzEz (Q/47reo) of the electric field, the induced 
super-spin-current Js (/ìq^Q/Stt), and the equivalent charge 
Qi (i.e. dzEl) versus the radiai distance r. 

the z-direction, then (s • \/)E^ = d^E^ = (^f+~2)5/2 
with = x"^ -\- y'^. Solving Eq.©, one has the induced 
super-spin-current density Js = — ^4^^ (^^2^^2)3/2 • This 
Js flows along the tangential direction (see Figl2fa)) and 
its spin points to the z-direction. On the other hand, 
as a usuai spin current,[20] the super-spin-current den- 
sity Js can generate an electric field E* in space which 
is the same as that generated by the electric dipole mo- 
ment Pe OC (— (^2^^2)3/2 7 Q7 Q) oi" the equivalent charge 

Qi = — V-pe c<: ^^2^^)5/2 • III FigISfb), we plot the radiai 
distributions in graphene for the variation dzE^ of the 
electric field of the originai charge Q, the induced super 
spin current density Jg, and the equivalent charge Qi (i.e. 
the spatial variation dzEl). For r < V^Z (r > V2Z) 
with dzEz being negative (positive), the spatial variation 
dzEl of the electric field E* induced by the super-spin- 
current is positive (negative). As a result, dzEl counter- 
acts the variation dzEz and then causes the variation of 
the total electric field in the SSC to vanish. 

Josephson effect is another highlight of superconduc- 
tivity and has wide applications.[3] We now investigate 
the similar effect for this SSC by considering a device con- 
sisting of two SSCs which are weakly coupled by a normal 
conductor, i.e., a SSC/normal conductor/SSC junction 
(see Fig. [sfa)). We can explicitly show the existence of 
the super-spin-current in this device in equilibrium. The 
weakly coupled junction is described by the Hamiltonian 
H = E/3(/3=L,i?) Hp^Hc^ Ht, where 




He = ^{ed^aMdHc^, 

a 



Namely, H^/r^ He-, and Ht are the Hamiltonians of the 
left/right SSC, the normal conductor, and the tunnel- 
ings between them, respectively. The order parameters 
^L/i?k = Ae*"^^/^, where the SSC phases (j)L/R are as- 
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FIG. 3: (color online), (a) The schematic diagram of the 
weakly coupled SSC-normal conductor-SSC junction. (b) The 
spin current Is vs. the phase difference A(j) for different A 
and Md = 0. (c) The Is-^c/) curves for different rud and 
A O.IM. Here = and T = 27rt|p/e O.IM; pk is the 
density of state of FM graphene in moment um space. 



sumed to be independent of the momentum k. We con- 
sider a phase difference Acj) = — (t^R between the left 
and right SSCs, which origins from a spin current flow- 
ing through the junction under the drive of an external 
device or from a variation of an external electric fìeld 
thread the ring junction device. The normal conductor 
is described by a level (or a quantum dot) with the spin 
index cj and spin-split energy Md- 

The spin-dependent particle current I^cr with the spin 
cr from the j3 SSC to the centrai normal conductor 
can be calculated by the following equation.|2l| /^^r = 
Re{2t*^/h) J ^G^cr ai^)^ where the lesser Green function 
G^^^(e) is the Fourier transformation of G^^ ^{t) = 

i(cj.(0)a/3k±cr(^))- This Green function G^^^(e) can be 
calculated by using the Dyson equation and so is the 
particle current //3cr.[21] Therefore, one can obtain the 
spin current Is = {Ir^ — lR^)h/2 and the charge current 
le = (^i?t + ^Ri)^' The charge current /g is identically 
zero because the e-h pairs are charge neutral. The spin 
current Is versus Acj) in the equilibrium with zero bias 
and zero spin bias is calculated and shown in Fig. [3l 
There is a super-spin-current flowing through the junc- 
tion that resembles the Josephson tunneling in a conven- 
tional superconductor junction. While Fig. [3l^b) exhibits 
the /s-A(/) curves for different values of the gap A with 
rrid = 0, Fig. [3fc) shows Is can also be observed in non- 
zero Md as long as A / and A^ ^ 0, tt. 

In conclusion, a SSC state as well as the electric 'Meiss- 
ner effect' and spin-current Josephson effect are predicted 
in FM graphene. For detection of the SSC, one can mea- 
sure the zero spin resistance or super spin current. Based 
on the experiment in Ref.[12|, here we propose a four- 
terminal device (as shown in Fig. le) which can be used 
to measure the non-local resistance and then confìrm the 



SSC state. [22^ A non-local resistance measurement means 
that a current is applied to the two right electrodes and 
the bias is measured on the two left electrodes. [12| In 
this device, the non-local resistance is induced by a pure 
spin current flowing from right side through FM graphene 
(SSC) to left side.[12| Once FM graphene turns into the 
SSC state by lowering temperature T < Te, the spin cur- 
rent will flow with no dissipation, leading to a sharp in- 
crease of the non-local resistance (see the Fig. le in the 
supplemental material). Due to the zero spin resistance 
effect in SSC, this sharp increase can be observed even 
in a macroscopic SSC device. 
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Supplementary information: 

In this supplementary material we discuss how to de- 
tect the spin superconductor state. In the following we 
first suggest fìve measurable physical quantities or meth- 
ods, and then follow up by proposing an experiment al 
setup. 

1) When the system enters the spin superconductor 
state, an energy gap opens up (see Fig.ld in the paper). 
This energy gap can be measured by ARPES or STM. 
When the temperature T is lower (or higher) than the 
criticai temperature Te, the gap emerges (or disappears). 

2) Because of the opening of an energy gap and the 
spin superconductor state does not carry charge current, 
the resistance sharply increases when the FM graphene 
enters from the normal metal state to the spin supercon- 
ductor state. This sharp increase in resistance can be 
experimentally tested. 

Not ice that the measurement s 1) and 2) are commonly 
done. These measurements demonstrate the opening of 
an gap and give strong hint that the sample may enter 
into the spin superconductor. 

3) The zero spin resistance is a main characteristic for 
the spin superconductor. Due to the zero spin resistance, 
the spin current can flow without any dissipation even for 
a macroscopic sample. At the end of the supplementary 
material we propose a four-terminal device and in this 
device one can show zero spin resistance by measuring a 
non-local resistance. 

4) In the Meissner effect for spin superconductor, an 
electric fìeld applied to the spin superconductor can in- 
duce the super spin current on the surface of the sample. 
This induced super spin current can generate an electric 
field E* that against the variation of the external electric 
field. Here the induced electric field E* is equivalent to 
that generated by a certain surface charge distribution. 
In this case, the induced electric fìeld and the equivalent 
surface charge distribution are the measurable quantities. 
For the example of Fig. 2 in the paper, the equivalent sur- 
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face charge Qi is proportional to ^ (^^2^^2)5/2 7 this value 
is quite large due to the smallness of the effective mass m* 
in graphene. But the equivalent surface charge Qi is stih 
finite even if m* = 0, since once it reaches the complete 
shielding and the variation of the total electric field is 
zero, no more super spin current and charge are induced. 
By considering the complete shielding case, the inducing 
surface charge Qi = j^r^^^^^^^-, where d ^ O.lnm is 
the thickness of graphene. Let us estimate the value of 
Qi. We assume that the distance Z between the external 
charge Q and the piane of graphene is lOnm and Q is 
a basic charge. The inducing surface charge Qi is about 
lO^^m"^ at the position r = 0. This surface charge is 
quite large and should be measurable. 

5) The super spin current in the the spin super con- 
ductor state or in the spin-current Josephson effect is 
also a measurable physical quantity. For example, it can 
be directly measured by observing the second-harmonic 
generation of the Faraday rotation as done in reference. 
Nature physics 6, 875 (2010). 

In addition, many indirect methods have successfully 
measured the spin current, e.g., by probing the spin ac- 
cumulation (Science 306, 1910(2004)), or probing the 
bias due to the inverse spin Hall effect (Nature 442, 
176(2006)), or probing the bias through a quantum point 
contact (Nature 458, 868(2009), etc. These methods can 
also be used in our system. Let us imagine a spin cur- 
rent fiowing into the spin superconductor (FM graphene) 
from one terminal, this spin current can fiow through the 
spin superconductor with no dissipation and fiows out to 
another terminal. Then we can measure the (usuai) spin 
current at the outside through the aforementioned meth- 
ods. 

Notice that it is not necessary to do ali of the afore- 
mentioned measurements. In fact, if one can take a mea- 
surement either in 3), 4), or 5), it establishes the presence 
of the spin superconductor. 

Now let US suggest an experiment setup and 
propose a detailed experimental process. 

Based on the device in the paper of Nature 488, 571 
(2007), we propose a four-terminal device consisting of 
a graphene ribbon coupled by four FM electrodes and a 
long FM strip, as shown in the Fig.l(a) and (b). Ex- 
cept for the long FM strip (the red region), the device 
is what used in the paper of Nature 488, 571 (2007). 
So this device can be realized by the present technology. 
The ideal case is that the long red strip is a FM insula- 
tor which directly couples to the graphene without the 
AI2O2, layer. But even if this is hard to achieve, it is also 
fine for the red region being a FM metal. Graphene in 
the red region has a FM exchange splitting, and it is nor- 
mal FM graphene at high temperature and turns to the 
spin superconductor at low temperature. Now we can 
qualitatively analyze the measurement results if the spin 
superconductor is realized. To simplify the analysis, we 
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FIG. 4: (color online), (a) and (b) are schematic diagrams for 
the proposed device, (a) is the top view and (b) is the side 
view. (c) and (d) schematically show the non-local resistance 
and the resistance i?23 versus the temperature, respectively. 

only consider the magnetic moments in ali FM electrodes 
and strip are in the same direction. 

Here we measure the non-local resistance as in the 
experimental paper of Nature 488, 571 (2007). In this 
measurement, a current is applied to two right electrodes 
(electrodes 3 and 4) and the bias is measured on two left 
electrodes (electrodes 1 and 2). When a current is ap- 
plied between electrodes 3 and 4, it injects a pure spin 
current into graphene. This pure spin current fiows from 
the electrode 3 through FM graphene (spin superconduc- 
tor) to the left side (as shown in the Fig.l(b)). Then it 
induces the bias between the electrodes 1 and 2. When 
FM graphene is in the normal FM metal phase, it has a 
finite spin resistance. In this case, the spin current grad- 
ually decays along its transport direction, so the bias and 
the non-local resistance i^i2,34 are small. Denoting Lij 
the distances between the electrode i and j and assum- 
ing L12 and I/34 are much longer than the spin relax- 
ation length A of the graphene in the normal state, the 
non-local resistance i?i2,34 can be obtained analytically 
(see Nature 416, 713(2002); Nature 442, 176(2006); etc): 
^norrnai ^ _h_^-L^2> I \ ^ ^hcrc W is the width of the 
graphene ribbon, P the spin polarization of the FM elec- 
trode, and GQ the conductivity of normal FM graphene. 
On the other hand, when FM graphene is in the spin 
superconductor phase, the spin resistance is zero and 
the spin current can fiow through it with no dissipation. 
In this case, the non-local resistance can be obtained 

00 pSSC _ 1 p2 A -(L2:^-Lssn)/X _ pnormal s se / ^ 
^12,34 (tgW^ ~ ^12,34 ^ ' 

which are much larger than R^2^2a^K Therefore, a sharp 
increase in the non-local resistance can be observed (see 
the Fig.l(c)) when FM graphene turns from the normal 
metal phase to the spin superconductor phase as tem- 
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pelature is lowered. Furthermore, the spin resistance 
in the spin superconductor is always zero regardless of 
the length of the spin-superconductor device, so the non- 
local resistance is independent of the length of spin su- 
perconductor system and the sharp increase in non-local 
resistance can be observed even in a macroscopic spin 
superconductor device. Experimentahy, one can also ob- 
serve the relation between the non-local resistance and 
the sample length to confìrm the zero spin resistance. 

In addition, if the long red strip is a FM insulator, 
a sharp increase can also be observed in the resistance 
R23 between electrodes 2 and 3 (see the Fig.l(d)), as 
discussed in point 2). 

Finally, we emphasize that the main body of the pro- 
posed device as well the pure spin current injected into 
graphene have been realized in an experiment (see Na- 
ture 488, 571 (2007), etc). So the proposed method is 
experimentally feasible. 
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